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ABSTRACT 

We report the first detection of a time-dependent, dust-scattered X-ray halo around a gamma-ray 
burst. GRB 031203 was observed by XMM-Newton starting six hours after the burst. The halo 
appeared as concentric ring-like structures centered on the GRB location. The radii of these structures 
increased with time as t 1 / 2 , consistent with small-angle X-ray scattering caused by a large column of 
dust along the line of sight to a cosmologically distant GRB. The rings are due to dust concentrated 
in two distinct slabs in the Galaxy located at distances of 880 and 1390 pc, consistent with known 
Galactic features. The halo brightness implies an initial soft X-ray pulse consistent with the observed 
GRB. 

Subject headings: gamma rays: bursts — X-rays: general — Galaxy: structure — ISM:dust 



1. INTRODUCTION 

It has long been realized that the small-angle scatter- 
ing of X-rays by dust grains can result in a detectable 
X-ray "halo" around a distant X-ray source (Overbeck 
1965), with a radial intensity distribution that depends 
on the dust properties and location. This phenomenon 
was first detected by Rolf (1983) using data from Ein- 
stein and confirmed by later observations (e.g. Mauche 
& Gorenstein 1986; Predehl et al. 1991). 

In the case of Gamma- Ray Bursts (GRBs), the tran- 
sient nature of the burst combined with the large ini- 
tial X-ray flux means X-ray scattering may be visible 
for a short period of time with a relatively high surface 
brightness. As the transient flux passes through a dust 
slab between the observer and the GRB, the scattering 
process introduces a time-delay in which the X-rays at 
larger angles to the line of sight arrive increasingly de- 
layed with respect to the non-scattered X-rays. Thus, 
the dust allows us to view the GRB X-ray flux at earlier 
times. The time-dependent X-ray halo around the GRB 
can, in principle, be used to provide detailed informa- 
tion on the location, spatial distribution and properties 
of the dust and the distance to and brightness of the GRB 
(e.g. Trumper & Schofelder 1973; Klose 1994; Miralda- 
Escude 1999; Draine 2003). Images of transient X-ray 
sources seen in scattered light are analogous to "light 
echos" seen around some supernovae (e.g. SN 1987A, 
Xu, Crotts & Kunkel 1995). 

Scattered X-rays observed at an angle 6 from the line- 
of-sight to a GRB arrive after the direct emission with a 
time-delay r s given by: 

Electronic address: sav2@star.le.ac.uk 

1 Department of Physics & Astronomy, University of Leicester, 
Leicester LEI 7RH, United Kingdom. 

2 Laboratory for High Energy Astrophysics, Code 662, NASA 
Goddard Space Flight Center, Greenbelt Road, Greenbelt, MD 
20771, USA. 

3 Universities Space Research Association 

4 Niels Bohr Institute for Astronomy, Physics and Geophysics, 
University of Copenhagen, Julian-Maries Vej 30, DK-2100, Copen- 
hagen 0, Denmark. 

5 XMM-Newton SOC, Villafranca, Apartado 50727, E-28080 
Madrid, Spain. 



(l + z s )D s D g 6 2 



where z s and D s are the redshift and angular diam- 
eter distance of the scattering dust, D g is the angular 
diameter distance of the GRB and D gs is the angular 
diameter distance from the dust to the GRB (Miralda- 
Escude 1999). 

Here we discuss the case of GRB 031203 observed by 
XMM-Newton on two occasions shortly after the burst. 
The first observation revealed the first dust-scattered X- 
ray halo detected around a GRB. 

2. OBSERVATIONS AND DATA ANALYSIS 

GRB 031203 was detected by the IBIS instrument on 
Integral on 2003 December 3 at 22:01:28 UT as a single 
peaked burst with a duration of 30 s and peak flux of 
1.3 x 10~ 7 erg s" 1 cm" 2 in the 20-200 keV band (Gotz 
et al. 2003; Mereghetti & Gotz 2003). A 58 ks XMM- 
Newton observation of the field began at 2003 Decem- 
ber 4, 04:09:29 UT. Several X-ray sources were detected 
within the 2.5' radius Integral error-circle (Santos-Lleo 
& Calderon 2003), the brightest of which (SI) appeared 
to fade through the observation (Rodriguez-Pascual et 
al. 2003) and was interpreted as the X-ray afterglow of 
GRB 031203. We refined the astrometric position of 
the afterglow by performing a cross-correlation with the 
USNO-A2 catalog, the improved SI position was (J2000) 
RA=08 h 02 m 30.19 s , Dec= -39° 51' 04.0" (la error 0.7", 
Tedds et al. 2003). 

The position in the Galactic plane (/ = 255°, b — 
—4.6°) results in high optical extinction (E(B — V) = 1.0, 
Schlegel, Finkbeiner & Davis 1998). Optical/IR obser- 
vations failed to locate an afterglow. Radio observations 
did locate a transient source location consistent with the 
X-ray afterglow (Soderberg, Kulkarni & Frail 2003), and 
optical imaging and spectroscopy revealed a candidate 
host galaxy at z = 0.105 which shows star- forming fea- 
tures typical of GRB hosts (Prochaska et al. 2003a, b). 
Our refined location for the X-ray afterglow is within 
0.5" of the optical galaxy and consistent with the posi- 
tion of the radio source. 
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Fig. 1. — Combined MOS images of GRB 031203 covering the 0.7 — 2.5 keV range, spanning 10' on a side. The data were divided into 
10 contiguous time intervals lasting 5780 s. The images were obtained from the first six time slices are shown (smoothed using a 6 arcsec 
Gaussian kernel). The three brightest point sources (SI, S2, S3) are marked. 



During the first XMM-Newton observation the GRB 
was observed on axis with all EPIC instruments in full- 
frame modes (Striider et al. 2001; Turner et al. 2001). 
A second 54 ks XMM-Newton observation was obtained 
on 2003 December 6 in which the GRB was « 6' off 
axis. The observational details and the GRB spectrum 
and light-curve are discussed in detail in Watson et al. 
(2004). 

The extraction of science products followed standard 
procedures using the XMM-Newton Science Analysis 
System (SAS v5.4.1). The EPIC data were processed 
using the SAS chains to produce calibrated event lists 
and remove background flares. GRB afterglow source 
counts were extracted from circular regions (radius 34") 
and background counts were estimated from a large off- 
axis region free from obvious point sources and the halo. 
The GRB afterglow fades as (t — to) -0 ' 4 , where to is the 
time of the GRB, which is unusually slow for GRB X-ray 
afterglows. 

3. THE X-RAY HALO 

The X-ray image from the first XMM-Newton obser- 
vation reveals an extended, circular halo concentric with 
the GRB 031203 afterglow (see Fig. 1). This halo was 
seen in all three cameras of the EPIC instrument, is 
unique to this observation, and is not due to scattered 
optical or X-ray light within the instrument. Summing 
the entire observation the halo had the form of a virtually 
complete ring (Vaughan et al. 2003), but splitting the ob- 
servation into contiguous time intervals revealed distinct 
rings which increased in radius through the observation. 

In order to better quantify this expansion, EPIC MOS 
images were produced for ten contiguous time intervals 



of 5780 s duration, spanning the entire first observation, 
and binned to 4" pixels. We concentrate on the MOS 
data as the halo lies over several chip gaps in the EPIC pn 
images, but the same halo was seen in both. The radial 
profiles of the images about the afterglow were calculated 
removing the three brightest point sources within the 
region of interest (SI, S2 and S3 in Fig. 1) by ignoring all 
counts within 40" from the centroid of each source. The 
radial profiles were then background subtracted using the 
mean level from the 350 — 400" annulus, which is well 
outside the detected halo. The resulting count profiles as 
a function of radius are shown in Fig. 2. The strongest 
peak moved outward from ~ 120" to ~ 220" during the 
observation, while a second peak moved from ~ 160 to 
- 270". 

The radial positions of the two rings were measured 
from the local maxima of the two peaks in the radial dis- 
tributions. Fig. 3 shows the change in radii as a function 
of time, and both were well fit by a simple power-law 
(0 oc (t — to) a ). Fixing to = 0, the indices of the power- 
law expansions were measured to be a = 0.53 ± 0.04 and 
0.51 ± 0.05 for the inner and outer rings respectively (er- 
rors are 90% confidence intervals), consistent with the 
cx (t — to) 1 / 2 expansion predicted for a scattered halo. 
As a consistency check we allowed to to be a free pa- 
rameter but fixed a = 0.5 and refitted to derive limits 
of t = 2794±i£?| s and 2005^867 s for the inner and 
outer ring respectively, both consistent with to = 0. The 
second XMM-Newton observation obtained 3 days after 
the burst showed no evidence for any ringed structure. 

The distance to the GRB is large (adopting the red- 
shift of the putative host galaxy), the time-delay is short 
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Fig. 2. — Radial profile of counts about the GRB 031203 afterglow from ten contiguous time intervals (derived from the combined MOS 
0.7 — 2.5 keV images). The afterglow itself (SI) and the next two brightest points sources (S2, S3) were removed prior to calculating the 
profiles (the positions of S2 and S3 are marked in the upper panels). The two peaks corresponding to the two expanding rings are obvious. 
The dotted lines denote the best-fitting oc (t — to) 1 / 2 functions (see Fig. 3). 



and is relatively large, thus to a good approximation 
D g = D gs and (1 + z s ) = 1 in equation (1), hence 

D s = 2cr s /0 2 . Assuming oc (t - t ) 1/2 and t = 0, the 
data in Fig. 3 imply that the scattering dust slabs are 
located at distances from the observer of D\ — 1388 ±32 
pc and D2 = 882 ± 20 pc, corresponding to the inner and 
outer ring respectively. At these distances, the ring di- 
ameters probe size scales of 2-3 pc in the dust scattering 
medium. 

The halo spectrum was extracted from an annulus with 
radii 110 — 220" using only the first half of the obser- 
vation; the decreased surface brightness of the halo at 
later times makes spectral extraction unreliable. A back- 
ground spectrum was extracted from a large off-axis re- 
gion (avoiding the halo and other point sources). Fig. 4 
shows the halo spectrum compared to the afterglow spec- 
trum. The spectrum of the halo is much steeper than 
that of the afterglow. Fitting with a simple absorbed 
power-law yielded photon indices of T = 1.98 ± 0.05 
and T = 3.03 ± 0.14 for the afterglow and halo, re- 
spectively. The best-fitting absorbing column was TVh = 
8.8 ±0.5 x 10 21 cm -2 when fixed to be the same for both 
afterglow and halo spectra. This is larger than the Galac- 
tic value of TVh = 6.1 x 10 21 cm -2 (Dickey & Lockman 
1990). The specific fluxes at t = 36650 s after the burst 
at 1 keV are 1.5±o.i x 10 -4 photons cm -2 s _1 keV" 1 
for the afterglow and 6.0 ±0.4 x 10 -4 photons cm -2 s _1 
keV -1 for the halo. 

To model the dust scattering medium, we have simul- 
taneously fitted the ten radial profiles shown in Fig. 2 us- 
ing the differential and total scattering cross sections of 



the Rayleigh-Gans approximation (Mauche and Goren- 
stein 1986). This is a valid approximation if the grain 
radius in fim is much less than the X-ray photon energy 
in keV (Smith & Dwek 1998). We assume that the dust 
is confined to two slabs at distances D\ and D2 with 
thickness AD. The halo seen is created by the com- 
bination of the soft X-ray pulse from the GRB and the 
subsequent afterglow. The observed angular width of the 
rings results from a combination of the duration of the 
soft X-ray pulse, At, AD and the point spread function 
(PSF) of the telescope. The width of the PSF corre- 
sponds to At - 1000 s or AD - 100 pc. Given that the 
GRB is very short, - 30 s (Mereghetti & Gotz 2003), 
we assume that At <C 1000 s and hence, allowing for the 
PSF, the observed broadening of the rings, 20", gives us 
an estimate of the slab thickness, AD = 130 ± 50 pc. 

The angular and temporal distribution was well fitted 
using grains with radii in the range 0.15 < a < 0.25 
/im. Smaller grains scatter to much larger angles (not 
visible because the time elapsed since the burst is too 
small) and the upper limit to the grain size distribution is 
expected to be ~ 0.3 /am (Predehl et al. 1991). Using this 
range of grain sizes, we require the initial GRB pulse and 
afterglow to have a photon index of T ~ 2.0 (as observed 
for the afterglow from 6 h) to match the observed photon 
index of the halo (r = 3.0). 

The brightness of the halo depends on the product of 
the differential scattering cross sections of the grains, the 
column densities of the grains in the slabs and the inte- 
grated flux of the X-ray pulse and afterglow. Interstellar 
extinction maps show structure in the GRB direction, 
with a considerable increase to Ay ~ 2 mag at a dis- 
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Fig. 3. — Expansion of the two rings around GRB 031203 
with time. The radius of each ring was measured from the local 
maxima of the radial profiles (Fig. 2). In both cases the expansion 
was well- fitted with a functional form (x. (t — to) 1 / 2 . The errors 
were assumed to be the width of the bin size used for the radial 
profile measurement. 
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Fig. 4. — Absorbed power-law model fits to the EPIC pn spectra 
of the GRB 031203 afterglow (circles, thick line) integrated over 
the entire observation and the dust halo (squares, thin line) inte- 
grated over the first half (27895 s) of the observation (during which 
the surface brightness was high enough for a reasonable spectral 
extraction). This clearly demonstrates the spectrum of the halo is 
significantly steeper than the afterglow. 



tance of ~ 1.3 kpc (Neckel & Klare 1980). This distance 
is consistent with the location of the more distant dust 
slab we detect. Adopting a mean grain size of a = 0.2 
/im, an Ay = 2 mag corresponds to a grain column den- 
sity of 1.5 x 10 8 cm -2 (Mauche and Gorenstein 1986). 
The best fit to the ring structure (x 2 / d -°- f - = 836/489) 
results from dividing this column density between the 
two slabs, Ni ~ 1 x 10 8 and N 2 ~ 0.5 x 10 8 cm" 2 (i.e. 
most of the dust is in the more distant slab). Using these 
values and a = 0.20 ± 0.05, the time-integrated flux re- 
quired for the soft X-ray pulse is 1600 ±800 photon cm -2 
keV -1 at 1 keV. Assuming a 30 s rectangular pulse this 
is equivalent to ~ 5 Crabs; using T = 2 the predicted 
GRB flux in the 20 - 200 keV band is 2.4 ± 1.2 photon 
cm -2 s _1 . Given the assumptions (e.g. spectral shape) 
this is in reasonable agreement with the value reported 
for GRB 031203 (peak of - 1.2 photon cm" 2 s" 1 , Gotz 
et al. 2003; Mereghetti & Gotz 2003). The X-ray burst 
is very bright compared to the afterglow (see Watson et 
al. 2004). Extrapolating the decay of the afterglow (sec- 
tion 2) to the time of the GRB accounts for only ~ 2 per 
cent of the halo emission. 

The line of sight to GRB 031203 passes close to the cen- 
ter of the Gum Nebula which appears as a 28° diameter 
sphere in Ha centered at I — 258°, b = —2° (Chanot & 
Sivan 1983). The Gum Nebula is likely to be a supershell 
created by repeated SNe explosions within it. Reynoso & 
Dubner (1997) detected a neutral gas disk (H i) associ- 
ated with the Gum Nebula, located ~ 500 ± 100 pc from 
us, with a radius of ~ 150 pc. The dust slab nearest to 
us may be associated with the rear face of the Gum Neb- 
ula which is known to have a concentration of molecular 
clouds (Woermann, Gaylard & Otrupcek 2001). 
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